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摘 要 : 为 实现 无 人 装备 在 野外 环境 下 ， 对 非 结 构 化 道 了 自动 、 普 适 和 精准 的 识别 与 导向 ， 提 出 了 一 种 基于 图 推 
模型 与 智能 寻 优 的 野外 场景 道路 导向 算法 。 首 先 将 图 像 分 割 为 同 质 超 像 素 块 ， 对 超 像素 块 的 多 特征 进行 融合 ， 构 造 训 
练 集 ; 改进 传统 拉 普 拉 斯 支持 向 量 机 算法 ， 结 合 超 像素 块 位 置信 息 动 态 选 取道 路 区 域 超 像素 种 子 块 ， 刘 央 是 邓 采 遇 
多 类 别 分 类 回归 器 和 相 邻 超 像素 的 一 致 性 回归 器 ; 结合 两 种 回归 器 的 回归 值 构造 马尔 科 夫 随机 场 的 能 量 函 数 ， 再 利用 
标准 图 割 算 法 选 代 求 取 最 小 化 能 量 函 数 实现 初始 道路 推理 分 结合 道路 初 分 割 结果 , 依据 人 对 道路 的 直观 感知 设 定 
约束 条 件 构 造 目 标 函 数 ， 利 用 差分 免疫 克隆 进化 算法 智能 寻 优 提取 道路 的 导向 线 。 在 南京 珠 山 采 集 的 数据 和 DARPA 
Grand Challenge 数据 库 上 进行 检测 ， 并 与 经 典 算法 的 道路 导向 效果 进行 定性 和 定量 比较 ， 结 果 表 明 该 算法 在 野外 环境 
下 对 非 结 构 化 道路 的 导向 线 提取 精度 总 体 达 91.79% 以 上 ， 相 比 于 经 典 算法 ， 检 测 精准 度 分 别提 升 48.1% 和 35.5%， 算 
法 处 理 效 率 分 别提 升 98.6% 和 97.8%， 在 检测 的 实时 性 与 精准 度 问题 上 实现 了 平衡 ， 具 有 和 较 强 的 应 用 前 景 。 
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Abstract: In order to realize automatic, universal and accurate identification and guidance of unstructured roads for unmanned 
equipment in the field environment, propose a road guidance algorithm for field scenes based on graph reasoning model and 
intelligent optimization. Firstly, segmentting the image into homogeneous superpixel blocks, and fusing multi-features of the 
superpixel blocks to construct a training set. Improving the traditional laplace support vector machine algorithm, combining the 
location information of superpixel blocks, dynamically selecting superpixel seed blocks in road areas, and training multi-class 
classifier regressors of superpixel blocks and consistency regressors of adjacent superpixels; Combining the regression values 
of two kinds of regressors, constructing the energy function of Markov random field and then using the standard graph cutting 
algorithm to iteratively obtain the minimized energy function to realize the initial road reasoning segmentation. Combining the 
initial segmentation results of roads, the objective function is constructed according to the constraints set by people's intuitive 
perception of roads, and using the differential Immune clonal evolution algorithm to intelligently optimize and extract the guide 
lines of roads. The data collected in Zhushan, Nanjing and DARPA grand challenge database are tested, and comparing the 
results qualitatively and quantitatively with those of classical algorithms. The results show that the extraction accuracy of the 
guide line of unstructured roads by this algorithm in the field environment is over 91.79 %, compared with classical algorithms, 
the detection accuracy is increased by 48.1 % and 35.5 % respectively, and the processing efficiency of the algorithm is increased 
by 98.6 % and 97.8 % respectively, which balances the real-time performance and accuracy of detection and has a strong 
application prospect. 
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其 基本 思想 是 双 直 线 反 映 道 路 的 大 致 方向 ， 
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有 具有 重要 的 现实 意义 和 广阔 的 应 用 前 景 。 


型 四、 埋 


估计 消失 点 与 道路 边缘 来 提取 导向 线 。 采 用 


基于 机 器 学 习 [5< 


见 阶段 的 道路 检测 的 算法 主 


要 分 为 基于 特征 1、 
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也 有 道路 标志 线 的 结构 化 道路 环境 已 经 较 好 地 完成 了 导 

些 特殊 的 如 军 寻 
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而 且 对 阴影 、 


向 线 提取 精度 。 但 是 对 于 边缘 纹 到 
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F 模 


三 大 类 。 文 献 [7] 提出 通过 迭代 的 方式 


偏差 较 小 。 
1 ”起 像素 类 标 MRF 推理 场景 分 割 


采 ) 


并 尽 可 能 大 地 涵盖 
或 ， 符 合 人 眼 对 道路 的 直观 感知 ， 可 操作 性 强 ， 出 现 的 


无 监督 超 像素 MRF 推理 类 标 技 术 进 行 野外 复杂 环境 


下 的 道路 初 分 割 。 既 考虑 到 超 像 素 块 的 多 特征 融合 信息 ， 而 且 


兼顾 了 超 像 素 块 领域 类 标 一 致 性 的 先 验 知识 
素 块 被 道路 内 超 像素 块 包 国 
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务 网 络 (vanishing point guided network,VPGN )， 联 合 处 理 
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且 在 非 结构 化 道路 条 件 下 效果 
针对 以 上 算法 的 不 足 ， 提 出 了 一 种 基于 
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| 二 倍 角 正弦 函数 对 方向 一 致 率 (orientation consistency 


里 方 


循环 过 代 的 方式 更 新 消失 点 和 边缘 ， 提 
高 了 消失 点 估计 精度 和 边缘 检测 精度 。 文 献 [8] 提 出 采用 纹 
单一 特征 的 局 


限 ， 


消 


和 识别 ， 具 有 较 强 的 鲁 棒 性 


在 


分 割 和 进化 算法 寻 优 的 野外 复杂 道路 导向 技术 。 
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像素 快 的 多 特征 融合 信息 ， 又 六 


的 先 验 知识 ， 若 某 些 疑似 超 像素 块 被 道路 内 超 像素 块 包 上 


技术 进行 野外 复杂 环境 下 的 道路 初 分 割 , 既 考 虑 到 超 
顾 了 超 像素 块 邻 域 类 标 一 致 性 
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加 倾向 于 分 类 成 道路 类 ; 若是 与 道路 整体 区 域 分 离 ， 更 加 倾向 


于 分 类 为 非 道 路 类 ， 有 助 于 提高 算法 


的 抗 


扰 能 力 和 重 棒 性 


更 加 符合 工程 实际 情况 。 第 三 ， 在 道路 类 聚 类 中 心 ( 超 像素 种 
子 块 ) 的 选取 上 融合 道路 模型 与 车 育 


二 沁 旨 


所 内 ， 必 定 为 路 的 先 验 知识 (参考 1 


1 最 小 左右 转弯 半径 交叉 履 
人 下文 ) 动态 选取 道路 类 


超 像素 种 子 块 , 不 仅 提升 了 算法 对 野外 “车 走路 变 ” 的 普 适 性 ， 


立 的 树叶 、 阴 影 、 
加 符合 工程 实际 情况 。 


， 如 某 些 疑似 超 像 


， 更 加 倾向 于 分 类 成 道路 类 ， 反 之 
助 于 降低 道路 区 域 阪 
水 迹 等 干扰 因素 对 分 类 结果 产生 的 影响 ， 更 


像 划分 为 同 质 


1.1 超 像素 块 多 特征 融合 
超 像素 将 图 像 从 像素 级 划分 成 区 域 级 ， 将 图 
的 区 域 再 进行 分 类 可 以 提高 图 像 分 割 的 效率 。 采 | 


i] SLIC(simple 


linear iterative cluster) 算法 进行 超 像 素 块 分 制 。SLIC 算法 与 传 


统 超 像 素 分 割 方法 相 比 处 理 速度 更 快 ， 


5 用 内 存 更 小 ， 边 缘 吻 


合 度 更 高 ， 可 以 将 图 像 划分 为 均匀 的 小 块 区 
容易 表达 ， 且 保留 物体 的 边缘 和 轮廓 等 重要 

为 了 获得 区 分 性 较 好 的 视觉 特征 ， 
述 超 像素 块 。 结 合 野外 非 结构 
文 提 取 以 下 四 类 特征 构成 视觉 特征 集 。 

a) 颜色 特征 。 野 外 道路 场景 图 像 的 颜色 
般 天 空 会 呈现 天 蓝 色 、 灰 白色 或 者 
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通过 多 种 特征 融合 的 方 
化 道路 场景 图 像 特点 ， 本 
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色 ， 道 路 则 是 褐色 或 者 红 
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褐色 。 只 有 寻 
色彩 模型 是 由 亮度 L 和 有 关 色 彩 的 a、 
于 感知 均匀 性 ， 它 的 工 分 量 密 
和 光照 有 较 强 的 鲁 棒 性 。 所 以 本 文 在 HSV ( 
度 ) 和 Lab 这 两 种 颜 是 取 颜 
在 Lab 颜色 空间 下 
(variance )、 供 
直方 图 和 饱和 度 

b) 纹 理 特 征 。 野 乡 


色 空 间 下 
个 颜色 


直方 图 四。 
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取得 很 好 的 区 分 效果 。 
LBP) 对 图 像 局 部 纹理 特征 具有 卓越 


还 有 效 避 免 了 随机 选取 道路 类 超 像素 种 子 块 可 能 引起 的 训练 仿 


差 和 效率 降低 。 第 四 ， 将 统计 学 习 理 论 与 最 优化 方法 


习 结合 起 来 ， 提 高 识别 算法 的 精度 条 
支持 向 量 机 (Lap SVMD) 的 


拉 斯 


规模 图 像 分 类 的 


I 运行 效率 。 针 对 
监督 分 类 方法 不 能 有 效 处 理 大 
问题 ， 提 出 了 通过 融合 道路 模型 与 4 
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基于 拉 普 


E 验 知识 动 
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空间 复杂 度 。 最 后 ， 借 助 统计 学 习 到 


有 效 地 提高 了 分 类 器 的 分 类 精 
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度 ， 降 低 了 算法 的 时 间 和 


和 最 优化 方法 解决 机 器 


通道 a、b 的 均值 (mean)、 
度 (skewness), 以 及 在 HSV 颜色 空间 下 


直 物 的 颜色 变化 比较 大 , 但 是 绿色 最 为 常见 。Lab 
b 三 个 要 素 组 成 ， 致 力 
团 匹 配 人 类 对 亮度 感知 ， 对 阴影 


色调 ， 饱 和 度 ， 明 


色 统 计 特征 。 提 取 超 像素 
方差 
的 色 度 


非 结 构 化 道路 场景 的 纹理 信息 体现 出 多 、 
杂 、 乱 的 特点 , 传统 统计 、 几 何 、 模型 等 方法 提 


取 纹 理 特征 不 能 


局 部 二 进 制 模 式 (local binary patterns， 
生 绘 能 力 , 计算 复杂 度 低 ， 


灰 度 尺度 不 变性 ， 易 于 工程 实现 ， 纹 元 字典 
以 非常 灵活 地 适应 计算 机 视觉 领域 。 其 中 


无 须 训练 学 习 ， 可 
MRELBP(median 


robust extended local binary pattern) 方 法 可 以 


有 获 微观 结构 和 宏 


观 织 构 信息 ， 具 有 计算 复杂 度 低 、 特 征 引 


数 较 低 、 对 高 


斯 随机 


AS 


噪声 、 随 机 像素 损毁 和 图 像 模 糙 


高 鲁 棒 性 
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噪声 、 椒 盐 
[93。 采 用 MRELBP 法 来 提取 超 像素 块 的 纹 


联 LBP 特征 谱 统 计 直 方 


里 特征 ， 


图 作为 特征 向 量 用 于 分 类 识别 。 


将 获得 的 关 


c) 结构 信息 。 结 构 信息 具 


有 很 好 的 稳定 


性 , 不 容易 受 外 界 
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录用 稿 华 夏 ， 等 : 基于 图 推 模型 与 智能 寻 优 的 野外 道路 导向 技术 


环境 的 和 干扰， 能够 有 效 提 高 超 像素 块 分 类 的 鲁 棒 性 和 抗 干扰 能 ” 图 像 中 尺寸 ，R 为 汽车 最 大 转弯 半径 对 应 图 像 中 尺寸 。 
力 。 采 用 稠密 方式 提取 尺度 不 变 特 征 变换 (scale-invariant feature 1.3 Lap SVM 超 像素 块 回归 器 训练 
transform, SIFT〉 和 方向 梯度 直方 图 (histogram of oriented 
gradient, HOG )。 即 先 在 超 像素 内 计算 每 个 像素 的 SIFT、HOG 
特征 向 量 ， 再 计算 超 像素 块 内 所 有 SIFT、HOG 向 量 的 几何 平均 
作为 其 特征 向 量 吕 。 
d) 位 置 和 形状 特征 。 野 外 非 结 构 化 道路 场景 中 超 像 素 块 的 
坐标 和 形状 能 够 为 超 像素 块 分 类 提供 很 强 的 几何 布局 线索 。 


超 像素 w 的 全 部 特征 向 量 记 为 有 ， 类 标 集合 表 示 为 


1 € {1,2,3} ， 分 别 对 应 道路 、 垂 直 障 碍 和 天 空 三 种 类 别 。 由 


广 和 1 组 成 多 类 别 回归 器 玖 ,的 训练 集 $ 如 式 〈2) 所 示 。 


9 ={ 太 0 O) 

超 像素 块 的 领域 信息 能 够 辅助 对 该 超 像素 的 属性 进行 判别 ， 
因此 训练 了 回归 器 用 来 估计 相 邻 超 像素 类 别 的 一 致 性 0。 一 致 
性 回归 器 采用 两 个 超 像素 的 差 值 进行 训练 09。 超 像素 对 的 一 臻 


~ 


图 1 SLIC 超 像素 图 像 分 块 

结果 可 以 看 出 ， 因 为 道路 区 域 纹理 、 颜 色 等 多 种 
杂 ， 所 以 分 割 形成 的 超 像素 块 形状 并 不 均匀 ， 与 天 
空 、 水 坑 等 特征 均匀 区 域 的 均匀 六 边 形 超 像素 块 对 比 明 显 。 采 
用 基于 Hu 不 变 矩 03] 的 方法 提取 形状 特征 , 即 归 一 化 的 超 像素 。 式 (3) 所 示 。 
不 变 矩 和 离心 率 值 。 
依据 一 定 的 模型 和 先 验 知识 ,针对 非 结构 化 道路 识别 问题 ， 
车 辆 在 行驶 过 程 中 ， 在 通道 宽度 内 ， 当 视觉 系统 安装 固定 且 沿 
车 辆 轴 向 时 ， 人 恒定 视野 区 域 始 终 为 可 靠 的 道路 区 域 09。 所 以 超 


性 标签 表示 为 有 e{1,1} ,这 里 的 -1 和 1 分 别 表示 相 邻 两 个 超 


像素 不 同和 相同 的 分 类 结果 ,一 致 性 


五 


归 器 瓦 , 的 训练 集 9。 如 


5 Fk) @) 


其 中 :8( 太 , ) 表示 相 邻 两 个 超 像素 的 特征 差 。 


像素 块 的 位 置信 息 同样 具有 较 好 的 区 分 度 ， 采 用 归 一 化 的 超 像 拉 普 拉 斯 支持 向 量 机 (LapSVM) 算 法 是 一 种 基于 流 形 正则 
素 块 中 心 像 素 位置 的 方法 提取 超 像素 块 位 置 特征 05。 的 半 监 督 分 类 算法 ， 研 究 如 何 同时 利用 少量 的 有 标志 样本 和 大 
1.2 道路 类 超 像素 种 子 块 动态 选取 量 的 未 标志 样本 进行 训练 和 分 类 P11。 通 过 引入 样本 流 形 正则 项 

针对 非 结 构 化 道路 识别 问题 ， 车 辆 在 行驶 过 程 中 ， 在 通道 20，LapSVM 算法 的 学 习 模型 中 将 包含 样本 的 固有 几何 结构 信 


宽度 内 ， 当 视觉 系统 安装 固定 且 沿 车 辆 轴 向 时 ， 人 恒定 视野 区 域 ” 息 。LapSVM 加 入 了 未 标志 样本 的 流 形 结构 信息 ， 并 将 这 一 信 
始终 为 可 靠 的 道路 区 域 。 在 道路 类 聚 类 中 心 〈 超 像素 种 子 块 ) 息 以 Laplacian 流 形 正则 项 的 方式 加 入 SVM。 学 习 过 程 中 充分 
的 选取 上 融合 了 车 前 最 小 左右 转弯 半径 交叉 覆盖 范围 内 ， 必 定 。 考虑 了 样本 间 的 局 部 几何 结构 ， 使 蕴涵 在 样本 中 局 部 的 鉴别 信 
为 路 的 先 验 知识 (参考 下 文 ;， 既 提升 了 对 野外 “车 走路 变 ” 的 。” 息 得 以 体现 ， 可 以 得 到 更 加 准确 的 分 类 结果 。 

自 适应 性 ， 同 时 又 避免 了 随机 选取 道路 类 超 像素 种 子 块 可 能 引 
起 的 训练 偏差 和 效率 降低 。 图 像 道 路 区 域 示 意图 如 图 2 所 示 。 


样本 集合 由 es {x,,1 =1,.…,n}, xX, 代表 第 i 个 样本 ,n 表 


示 样本 的 个 数 。 集 合 L={x,i=1…,m] 表示 标志 样本 ，m 表示 
标志 样本 的 个 数 ,U ={x,i=1,…,u} 表示 无 标志 样本 的 集合 ,4 


表示 未 标志 样本 的 个 数 。 用 y; 表示 第 i 个 样本 所 属 类 别 ， 


y E{ 一 ],1} 。LapSVM 模型 定义 为 


f =minZmax(l -yf 05)0) + | a 


图 2 图 像 道路 区 域 示 意图 

可 以 确定 道路 图 像 中 , 车 前 以 略 大 于 车 宽 的 尺寸 万 为 底 边 
长 , 高 为 1 的 等 腰 三 角形 区 域 , 为 图 像 中 道路 类 超 像素 种 子 块 选 
取 最 佳 区 域 。 依 据 几 何 学 知识 可 以 求 得 高 ! 如 式 (1) 所 示 。 


1= JR -(r+ZD2Z) (1) 


其 中 :万 为 车 宽 对 应 图 像 中 尺寸 ;7 为 汽车 最 小 转弯 半径 对 应 


| 


其 中 : 有 =[f (Xx;)，X; ES] ,在 训练 数据 集合 上 应 为 n 维 


列 向 量 。 | || 是 在 再 生 核 希 尔 伯 特 空间 (RKHS) 中 定义 的 环境 


范 数 (Ambient Norm)， 万 , 是 与 核 函数 相关 的 RKHS 。 参 数 y， 
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录用 稿 


是 权 


中 


E 来 控制 | 了, 在 RKHS 中 的 复杂 性 ，| 了 ,是 流 形 


j 来 保持 样本 分 布 的 内 在 流 形 结构 ， 


正则 项 (intrinsic regularizer), | 


参数 y) 是 低 维 流 形 中 函数 的 权重 , 来 控制 内 在 几何 结 


和 函数 的 


性 。 通 过 计算 拉 格 朗 日 乘 子 ， 得 到 分 类 器 为 


三 PoiK(, G) 


其 中 : w' 是 拉 格 朗 日 乘 子 ， 玉 是 核算 阵 。 求 解 拉 格 朗 日 乘 子 ; 


a =(277+27 KZ T° ,YB (0) 


其 中 :J 了 是 单位 矩阵 ; 


工 是 Laplacian 和 矩阵; 了 e R”™”" 是 由 标 


志 样本 y ,i=b...，m ， 组 成 的 对 角 和 矩阵 。J， EeE R”” 是 由 标 


志 样本 与 未 标志 样本 形成 的 分 块 矩阵 [L 0]。 p' 是 Lagrange 乘 


P=max DB -3p"0p 0) 
QO=Y,K(2yT+27y,KL) TTY (8) 


st.>By=0;0<pB sbi=b,m OO) 
i=] 


从 上 面 的 定义 可 以 看 出 , LapSVM 算法 中 涉及 许多 矩阵 的 运 
算 和 转换 ， 当 无 标记 样本 很 多 时 ， 需 要 很 大 的 内 存 空 间 和 很 长 
的 CPU 占用 时 间 , 甚至 引起 内 存 溢出 等 问题 。 可 以 通过 原始 优 


化 来 加 速 训练 过 程 ， 将 LapSVM 的 模型 式 重 新 定义 为 


QeR" ,beR 


m max(l—y,(K/g+b),0)+ya0 Ka 5 
Imin (10) 
+ (a K+I'bL(Ka+B)) 


其 中 : 了 为 SVM 定义 中 的 闵 值 ，y, ,i= 上 ...，m 为 标志 样本 ; 参 
数 y。， 是 权重 ， 是 核算 阵 ，w 是 拉 格 朗 日 乘 子 。 


利用 LapSVM 算法 进行 分 类 ， 标志 样本 数量 越 少 ， 则 分 类 
结果 误差 可 能 越 大 。 本 文通 过 上 文 动态 选取 道路 区 域 标志 样本 ， 
选取 那些 在 确定 道路 区 域 里 面 或 者 附近 的 超 像素 块 样本 即 可 能 
届 于 道路 区 域 的 超 像 素 块 样本 加 入 训练 样本 集 ， 舍 弃 了 离 确定 
道路 区 域 较 远 或 在 图 像 四 角 位 置 的 无 标记 样本 ， 使 得 参与 训练 
的 无 标记 样本 提供 更 丰富 的 启发 信息 。 

结合 动态 道路 预选 取样 本 的 Lap SVM 超 像素 块 分 类 算法 ， 
首先 对 全 部 无 标记 样本 选取 少量 含 启 发 信息 更 为 丰富 的 m 个 
无 标记 样本 。 假 设 有 /个 标记 样本 以 及 预选 取 m 个 无 标记 样本 ， 
以 1+ m 个 样本 建立 丰 近邻 图 然后 ， 在 此 图 上 构造 目标 函数 / 
来 实现 分 类 。 

算法 基本 流程 如 下 : 

a) 输入 标志 样本 集合 L={x,i=1…,m} 和 无 标志 样本 的 集 


合 U={x,i=1,…,U)}; 


华 夏 ,， 等 : 基于 图 推 模型 与 第 


b) 从 u 个 无 标记 样本 中 预选 取 含 启发 信息 更 为 丰富 的 m 
个 无 标记 样本 ; 
c) 利用 高 斯 核 函 数 ， 计 算 核 甜 K, = K(x,x) ; 


d) 计算 图 Laplacian 矩阵 L=D-W,D 是 对 角 和 矩阵 ， 


e) 选取 合适 的 权重 / ，y 


f) 计算 qr =(27 7+27 KDJT7D ; 


g) 输出 分 类 函数 广 = 立 w KG,x) 。 
i=l1 


通过 对 标志 样本 和 无 标志 样本 的 超 像素 块 特征 向 量 进行 训 
练 学 习 ， 可 以 得 到 Lap SVM 分 类 器 , 即 可 用 于 对 野外 复杂 非 结 
构 化 道路 图 像 进行 道路 区 域 检测 与 识别 。 通 过 训练 学 习 得 到 多 


类 别 回 归 器 及 。 在 测试 过 程 中 , 对 于 一 个 输入 特征 , 多 类 回归 


器 及 的 输出 是 1、2、3。 其 中 隐 含 层 节点 数 取 为 3， 可 以 自动 


地 动态 选取 道路 区 域 超 像 素 种 子 块 作为 道路 聚 类 中 心 。 


致 性 回归 器 已 2 的 回归 值 作 为 2 个 超 像素 块 属于 相同 


I 
运 


类 别 的 隶属 度 。 已。 回归 器 的 输出 也 是 -1 与 1 之 间 的 连续 值 。 


1.4 超 像素 类 标 MRF 推理 分 割 
马尔 可 夫 随 机 场 (Markov random fields) 模型 相 比 其 他 方 


法 的 优势 是 : 提供 了 一 种 简单 方法 来 对 先 验 知识 建 模 ,， 相 比 其 
他 基于 像素 点 ， 或 基于 局 部 的 方法 , 它 可 以 考虑 到 环境 知识 的 


五 


向 ， 如 果 建 立 的 图 模型 得 当 ， 进 而 可 能 获得 更 加 符合 人 眼 视 
局 最 优 解 09。 当 超 像素 分 类 完成 之 后 ， 采 用 MRF 推理 


误 法 
es 


的 方法 得 到 超 像素 的 类 标 。 对 于 一 幅 输 入 图 像 , 假设 A = {a} 


是 图 像 的 超 像素 集合 ,对 应 的 成 对 MRF 能 量 函数 如 式 〈11) 所 


No 


| 


~ 


EO=》 VC)+w > Vl) 


aieA (aisaj)eP 


(11) 


其 中 : i 和 j 是 超 像素 索引 号 ; [和 /分 别 是 超 像素 i 和 j 的 候 


选 标 签 ，Q 是 成 对 能 量 项 的 权 值 ，U'; 和 Vi 分 别 代表 数据 项 和 
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平滑 项 势能 函数 , 由 超 像素 的 


归结 果 计 算得 到 ，L/; 表示 分 配 


类 标 上 到 像素 块 a; 的 代价 值 ， 为 a; 的 类 别 回 归 值 与 类 标 值 之 


间 差 的 绝对 值 ， 如 式 (12)〉 所 示 。 


录用 稿 
U,(L)=|H(f) -1 (12) 


超 像素 对 {Qi,a)} 之 间 的 成 对 势能 V 是 惩罚 相 邻 超 像素 


4;、Q4; 取 不 同类 标的 代价 值 。 采 用 


Potts 模型 来 表示 成 对 的 


代价 函数 ， 如 式 (13) (14) 所 示 。 


V ;= BU IT,,L,) (13) 
,fz 
VOD | rar (14) 


其 中 : 8 是 相 邻 超 像 素 之 间 的 权 值 ， 
出 来 计算 ， 如 式 〈15) 所 示 。 


I 


类 别 一 致 性 回归 器 的 输 


Bs) =1+ H,(f,f)) C15) 


图 割 算法 (graph cuts) 是 优化 离散 方程 的 有 力 工具 ， 本 文 
考 文献 通过 最 小 化 图 割 来 最 小 化 MRF 能 量 函数 从 而 完成 对 
像 的 场景 分 割 [19201。 


E 寻 优 道路 导向 线 提取 


通过 道路 分 割 能 够 获得 复杂 不 规则 的 道路 区 域 ， 如 何 引导 
车 辆 下 一 步 的 前 进 方向 成 为 无 人 驾驶 的 又 一 关键 问题 。 对 不 规 
则 的 道路 边界 基于 纹理 消失 点 的 导向 线 估计 误差 较 大 ， 而 基于 
道路 模型 的 导向 线 估计 边界 在 复杂 野外 环境 下 待 拟 合 点 选择 困 
难 , 拟 合 偏差 较 大 。 双 直线 道路 边界 能 够 反映 道路 的 大 致 方向 ， 
并 尽 可 能 大 地 涵盖 道路 区 域 ， 可 操作 性 强 ， 出 现 的 偏差 较 小 。 
2.1 非 结构 化 道路 双 直 线 边 界 与 导向 线 模型 
大 部 分 野外 非 结 构 化 道路 并 不 是 直线 ， 但 通常 曲率 很 小 ， 
入 定 非 结 构 化 道路 的 左 、 右 两 边界 都 是 直线 ， 这 样 就 可 以 采 | 
直线 道路 边界 来 对 图 像 中 的 道路 边界 进行 约束 ， 如 图 3 所 示 。 
这 大 大 降低 了 运算 的 复杂 度 ， 提 高 了 程序 的 运行 效率 。 


Sp 


| 


2 ”智和 


图 3” 非 结构 化 道路 双 
依据 以 上 两 条 道 


线 边 界 及 导向 线 模型 
首 路 边缘 判定 准则 ， 本 文 以 图 像 的 左下 角 顶 


点 为 整个 坐标 系 的 原点 0， 直线 全 函数 方程 如 式 (16〉 作 为 道 


路 的 左边 缘 。 
=k xx+b (16) 


ChinaXiv 合 作 
等 : 基于 图 推 模型 与 箱 I 


二 
ba 


道路 的 右边 缘 取 为 直线 1/,， 函 数 方程 如 式 〈17) 所 示 。 


y, =k, xx+b, (17) 
片 的 尺寸 为 mkn， 由 方程 组 式 (18) 


pa 


»=kxx+b; 
y=k,xx+b,; 
0=k xx+D; 


18 
0=k, xx,+b; (18) 


可 求解 得 到 外 与 x 轴 交 点 C( Xi,0) ， 与 y 轴 交 点 B 


(07) 


”kk (19) 


(20) 
=b, 


二 与 1 交点 E( 义 ,也 ) 作为 非 结构 化 道路 的 消失 点 : 


CD 


依据 先 验 知识 假定 车 辆 当前 所 处 位 置 为 F 点 [2 ,0) ,由 


E、F 点 则 可 以 求 得 道路 的 导向 线 4 ,由 点 和 斜 公式 可 以 求 得 4 的 


方程 ， 如 式 (22) 所 示 。 
byek — bok 


x b -hk hk) 


n (bok — bek,) 
已 = 一 后. 
入 -> 2b, —2b, — ne(k —k,) 


» =kxxt+b 


22) 


车 辆 继续 前 进 需 调整 角度 为 Y， 如 式 (23〉 所 示 。 


tan(h) -Th ,hh<d) 
Th -tan(ks) ,hh>0 


(23) 


ChinaXiVE 3 
a 华夏 ， 等 : 基于 图 推 模 型 与 和 能 寻 优 的 时 外 道路 导向 技术 


改 、/( /三 十 日 工 1| 
一 人 下 上 日 二 
1, IF 


2.2 ”有 约束 道路 边界 双 直线 估计 目标 函数 构造 


非 结构 化 道路 的 双 直线 模型 结合 先 验 知识 ， 可 以 得 到 最 


结合 人 眼 对 道路 的 直观 感知 和 先 验 知识 ， 通 过 整体 估计 的 。 ”优化 目标 函数 的 约束 边界 条 件 : 

方式 来 进行 道路 边界 双 直 线 估计 。 a) 为 提高 运行 效率 ， 加 快 收敛 速度 ， 假 设 道路 消失 点 即 两 
合理 的 道路 边界 双 直 线 应 该 尽 可 能 完整 地 包含 道路 区 域 ，。 边界 交点 在 图 像 内 部 : 

即 边界 所 包含 的 道路 部 分 占 道路 总 区 域 的 比例 要 尽 可 能 大 。 同 a 

时 ， 最 优 的 道路 边界 包含 的 区 域内 ， 道 路 应 该 占有 尽 可 能 多 的 pr 

比例 。 由 此 构造 了 两 个 目标 子 函数 ， 边界 包含 区 域内 道路 部 分 Cy 

所 占 比例 L441、 边界 内 道路 占 道路 整体 的 比例 L 。 依 据 构造 的 

双 直 线 模型 可 以 得 到 双 直 线 模型 道路 边界 包含 的 区 域 面积 关 A 


SI: Be (10,170 )， 取 值 间隔 AQ = 5”， 则 两 直线 斜率 : 
Kk 10°), 170° 
bh bh i GD 
5 = (十 -一 )x( 一 +bl) /2 (24) k, < (tan(10°), tan(170°)) 
k h 4 h 


2.3” 群 智能 边界 寻 优 估计 
Dm 因为 道路 区 域 的 面积 一 般 是 不 规则 图 形 ， 通 过 初 分 割 可 以 差分 免疫 克隆 算法 (differential immune clone clustering 
之 获得 道路 边缘 像素 点 的 位 置信 息 ， 转 换 到 双 直 线 模型 坐标 系 中 algorithm,DICCA) 主要 是 基于 差分 进化 及 免疫 克隆 ， 通 过 克隆 
繁殖 、 差 分 变异 、 交 叉 及 克隆 选择 操作 来 进化 种 群 ， 在 进化 过 
程 中 加 入 局 部 搜索 机 制 来 提高 算法 的 收敛 性 。 DICCA 普 适 性 好 ， 


可 以 表示 为 坐标 (X;, y;), 使 用 牛顿 积分 求 边界 内 不 规则 道路 


区 域 的 面积 6。 不 易 陷入 局 部 极 值 , 有 较 好 的 全 局 收敛 性 和 鲁 棒 性 ,非常 适合 求 
解 各 种 数值 最 优化 问题 。 先 合理 构造 约束 条 件 和 目标 函数 ， | 

o=|do= Ddo ,Fig 约束 边界 条 件 和 最 优化 目标 函数 通过 DICCA 和 迭代 优化 的 方式 ， 

a 即 可 求 得 最 优化 道路 边界 的 参数 上 、b、k、b,。DICCA 算 


法 首先 需要 初始 化 参数 ,具体 有 种 群 规模 为 30、 交 叉 概率 为 0.85、 
(25) ”变异 概率 为 0.05、 最 大 迭代 次 数 为 100。 算 法 的 操作 算 子 有 差 
分 变异 、 差 分 交叉 、 克 隆 增殖 、 均 匀 变 异 、 克 隆 选择 [3 五 个 。 

经 过 克隆 选择 操作 ,得 到 下 一 代 种 群 。 算 法 的 终止 条 件 是 : 

在 种 群 欠 代 的 过 程 中 ， 当 达到 算法 最 大 和 迭代 次 数 ， 或 者 连续 4 
代 的 适应 度 值 都 小 于 最 小 迭代 的 精度 时 ， 算 法 停止 迭代 ， 输 出 

5 =》 | 一审 Xy，y=min(y,, yi) (26) 最终 隧 类 结果 ,算法 终止 .整个 DICCA 算法 具体 流程 如 图 4 所 
本 示 , 全 文 算法 流程 如 图 5 所 示 。 


Xx 


初 分 割 出 来 整体 的 道路 区 域 面积 5, : 


初始 化 种 群 2 
克隆 繁殖 操作 


均匀 变异 操作 


克隆 选择 操作 更 新 种 
群 2 


(ey 
5 (27) 
32 


变异 差分 操作 


本 文通 过 将 Wi 和 /2 线性 加 权 的 形式 构造 DICCA 算法 的 


标 函 数 F， 如 式 (28) 所 示 。 


1 


= 一 一 一 一 一 (28) 
NxXH + xH 


选择 操作 ， 更 新 种 群 1 


其 中 : 如。 、 克 分 别 代表 约束 条 件 子 函数 4 和 / 对 总 体 最 优 


用 种 群 2 中 的 个 体 取 
代 种 群 1 中 的 d 个 最 差 
个 体 


化 目标 函数 的 指标 权重 。 


所 求 得 的 最 优化 目标 函数 F 的 极 小 值 对 应 的 道路 边界 即 是 
最 优 检测 结果 。 


输出 最 优化 道路 边界 参数 


输出 新 种 群 


图 4 DICCA 迭代 寻 优 算法 流程 
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输入 野外 道路 场景 图 像 
SLIC 超 像素 分 块 


动态 选取 道路 超 像素 种 子 块 


超 像 素 块 多 特征 提取 ， 训 练 样本 构造 


RBF 训练 多 类 别 回归 器 二 ne RBF 训练 一 致 性 回归 器 Hl 


MRF 推 理 超 像素 类 标 图 像 分 害 


DICCA 优 化 导向 线 提取 


输出 野外 道路 场景 分 割 图 像 和 道路 导向 线 


图 5 本 文 算法 流程 


3 ”实验 结果 与 分 析 


为 了 验证 本 文 算法 对 野 儿 
测 的 有 效 性 和 导向 线 提取 的 准确 性 。 实 验 测试 了 多 种 典型 复 
非 结构 化 道路 ， 并 与 文献 [7] 中 的 算法 以 及 文献 [8] 和 文献 [9] 
的 算法 进行 识别 效果 比较 ， 比 较 算法 代码 均 来 源 于 网 络 开源 
源 。 对 实验 结果 采用 定性 和 定量 的 评价 方法 进行 分 析 。 
3.1 实验 设计 和 数据 库 


复杂 场景 下 非 结构 化 道路 边缘 检 


杂 


中 


次 
ba 


实验 以 四 核 Intel i7 处 理 器 ，8.0 GB 
台 , 在 Windows7 操作 系统 下 以 MATLAB R2014a 作为 开 
工具 。 实 验 所 用 非 结构 化 道路 图 像 数 
grand challenge ) 场景 分 割 数 据 库 和 在 南京 珠 山 采 集 的 非 结 构 
道路 数据 库 ， 以 及 互联 网 上 经 过 筛选 得 到 的 典型 复杂 非 结 构 
道路 图 像 ， 将 所 有 实验 图 像 进行 尺寸 归 一 化 为 640*480 像素 
3.2 ”实验 结果 评价 指标 

为 了 能 够 对 各 种 算法 的 道路 边缘 检测 质量 进行 定量 比较 
本 文采 用 文献 [7] 中 的 衡量 标准 ， 精 准 度 + 来 评价 各 算法 道路 


局 


内 存 PC 机 为 硬件 平 


发 


居 取 自 于 DGC (DARPA 


化 
化 


o 


9 


[x| 


域 分 割 精度 的 高 低 , + 越 小 检测 精度 越 低 ,反之 检测 精度 越 高 。 


BMB, 
ee 
BUB, 


其 中 万， 是 人 工 标定 的 道路 区 域 像素 集 ，B, 是 算法 检测 


的 道路 区 域 像素 集 ; 


局 B, 表 示 两 者 并 集 23。 


同时 ， 为 了 定量 评价 各 算法 的 导向 线 提取 精度 ， 本 文 定 
了 导向 精度 mm ， 如 式 〈33) 所 示 。 


1o0% -ek = maxesb G 


其 中 : k, 是 人 工 标注 的 导向 线 斜率 ; k, 是 算法 提取 出 的 导向 


(32) 


LH 
LI 


到 站 B, 表 示 两 者 交集 ， 即 共同 部 分 ; 


义 


3) 
线 


华 ， 夏 ， 等 : 基于 图 推 模型 与 智能 寻 优 的 野外 道路 导向 技术 


斜率 。 为 了 使 导向 精度 的 结果 更 加 符合 人 们 的 认 知 ， 对 z 进行 
数值 变换 转化 为 z,， 如 式 (34) 所 示 。 


=50%+ 辽 .100% (34) 


Z3 越 大 ， 导 向 线 提取 精度 越 高 。 


3.3 ”实验 结果 
指标 权重 的 合理 


性 直接 影响 着 多 属性 评价 结果 的 准 


为 了 合理 选 定 最 优 权 值 4 、 也 ， 提 高 道路 边缘 检测 精度 。 


为 了 降低 运算 复杂 度 ， 对 丸 、 刀 进行 数值 约束 。 


(35) 


在 数据 库 中 随机 抽取 20 幅 图 像 ， 人 工 标定 道路 区 域 和 道 


路 边界 直线 , 则 由 式 容易 求 得 J4 和 / 心 , 通 过 对 九 在 区 间 (0,1) 


以 间隔 0.001 取 值 循环 夫 代 ， 刀 =1 一 入 ， 使 得 目标 函数 F 取 


得 极 小 值 ， 记 录 下 对 应 的 4 、 见 值 进行 权 值 样本 的 构造 。 


03 
07 
05 
05 
~ 04 
03 
02 
01 
2 345678 9 1011121314151617 181920 
Image Number 
图 6 样本 数据 
通过 对 九 求 均值 的 方式 ， 获 得 了 普 适 性 较 好 的 最 优 权 值 


4 = 0.645 
] =0.355 


在 实验 中 定性 及 定量 地 对 比 了 本 文 和 文献 以 及 文献 提出 的 
非 结 构 化 道路 边缘 检测 算法 的 检测 性 能 ， 检 测 效果 如 图 7、 表 
1 一 3 所 示 。 
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2 可 以 发 现 ,在 道路 导向 精度 上 ， 本 文 算法 比 文献 [7，8] 分 别 高 
出 35.5% 和 48.1%， 比 文献 [10] 也 要 高 出 15.5%; 由 表 2 可 以 发 
现 , 在 算法 实时 性 上 ， 本 文 算法 比 文献 [7，8] 分 别 高 出 97.8% 和 
98.6%， 基 本 能 够 实现 道路 的 快速 处 理 。 

表 1 各 算法 道路 区 域 分 割 精度 


Tt 最 小 值 
算法 tT 均值 (100%) t 最 大 值 (100%) 
(100%) 
文献 7 61.86081 47.70483 96.48042 
文献 8 54.70859 35.50663 95.71827 
文献 9 91.76281 81.50393 99.96012 
本 文 算 法 89.89081 81.90483 99.68967 
表 2 各 算法 导向 精度 rs 
Tt3 最 小 值 
算法 Tt 3 均值 (100%) Tt3 最 大 值 (100%) 
(100%) 
文献 7 56.26182 0.51423 96.32143 
文献 8 43.70859 0.32643 85.41937 
文献 9 76.28165 69.59314 98.33963 
本 文 算法 91.79281 72.62583 99.79864 
表 3 各 算法 单 张 图 像 平均 处 理 时 间 
算法 文献 7 文献 8 文献 9 本 文 算 法 
运行 时 间 
28.32643 43.70859 0.0352 0.61937 
= 2 (S) 
图 7 各 算法 道路 分 割 与 导向 结果 
4 ”结束 语 


图 7 中 第 一 列 为 原 图 和 人 工 标注 道路 消失 点 ， 第 二 列 为 文 
献 [8] 算 法 检测 效果 ， 第 三 列 为 文献 [7] 算 法 检测 效果 图 ,第 四 列 本 文 提出 的 基于 图 推 模型 与 智能 寻 优 的 野外 场景 自 适应 道 
为 文献 [9] 算 法 检测 效果 ， 第 五 列 为 本 文 算法 检测 效果 。 由 于 文 ”路 导向 算法 。 实 验 结果 表明 本 文 算法 在 野外 复杂 环境 下 对 非 结 
献 [7] 和 文献 [8] 都 是 基于 像素 纹理 的 消失 点 检测 , 依赖 于 对 纹理 构 化 道路 的 道路 区 域 分 割 精度 总 体 达 到 89.9%， 导 向 精度 总 体 
村 征 的 提取 , 当 路 面 纹理 较 弱 不 清晰 如 图 6 中 第 4、10 张 图 像 ， 达到 了 91.8%, 算法 处 理 效率 分 别提 升 97.8% 和 98.6%, 综合 性 


或 者 杂乱 如 图 6 中 3、5、7、9 图 像 时 检测 效果 较 差 ; 文献 [10] 。 能 明显 优 于 经 典 算法 ， 在 检测 的 实时 性 与 精准 度 问 题 上 实现 了 

在 利用 非 结构 化 道路 训练 集 进行 训练 后 ， 对 光照 、 纹 理 等 干扰 ” 平衡 ， 具 有 较 高 的 道路 检测 精度 和 较 好 的 实时 性 ， 对 道路 区 域 

有 较 好 的 鲁 棒 性 ， 但 是 对 于 形状 不 规则 和 场景 混乱 的 非 结构 化 。 ”的 判定 和 道路 轮廓 走向 的 估计 更 加 符合 人 了 眼 对 道路 的 直观 感知 ， 

道路 识别 效果 较 差 ， 如 图 6 中 第 5、8、9 幅 图 像 ， 本 文 基于 经 。 ”能够 实现 对 野外 非 结构 化 道路 精准 快速 检测 ， 具 有 极 强 的 应 用 

过 初始 分 割 得 到 的 道路 区 域 再 进行 道路 边缘 检测 可 以 有 效 去 除 。 前 景 。 但 是 采用 的 双 直 线 道路 边界 估计 模型 较为 简单 对 于 弯 道 

无 关于 扰 ， 提 高 算法 的 运行 效率 ， 对 道路 边界 检测 的 效果 也 明 ”等 形状 复杂 道路 识别 效果 一 般 , 且 算 法 的 实时 性 仍然 有 待 提高 。 
bE 


i 


显要 优 于 其 他 两 种 算法 ， 能 够 有 效 克 服 强 光 、 阴 影 、 水 迹 、 雨 pe 

雪 、 雾 天 落叶 等 外 部 环境 对 道路 边缘 检测 的 干扰 ， 能 够 有 效 检 ， 

测 非 结构 化 道路 边缘 。 但 是 对 于 茶 些 形状 不 规则 ， 有 弯 道 的 非 。 [1] 王 晓 栋 , 徐 成 , 刘 彦 . 一 种 实时 鲁 棒 的 非 结构 化 道路 检测 算法 [用 . 计 

结构 化 道路 ， 本 文 算法 的 检测 仍然 存在 着 较 大 的 局 限 性 ， 如 图 算 机 应 用 研究 ,2010, 25 (12): 2763-2765. (Wang Xiaodong, Xu Cheng, Liu 

6 中 第 4、9、10、5 图 像 。 Yan. Real-time and robust method for unstructured roads detection [J]. 

为 了 降低 随机 性 和 偶然 性 对 分 割 性 能 评估 的 影响 ， 本 文 从 Application Research of Computers, 2010, 25 (7): 2763-2759. ) 

数据 库 和 互联 网 上 选取 了 总 共 100 张 非 结构 化 道路 图 像 ， 利 用 [2] 地 运河 ， 张 浩 峰 , 於 敏 杰 , 等 . 基于 K-means 特征 的 复杂 环境 下 道路 识 

各 算法 分 别 进行 道路 边缘 检测 和 导向 线 提取 ， 然 后 计算 道路 区 别 算法 [四 计算 机 应 用 研究 , 2016, 39 (2): 602-606. (Hao Yunhe, Zhang 
Haofeng, Yu Minjie, et al. Road recognition algorithm under complex 


~ 


域 分 割 精度 + 和 导向 精度 荆 ;， ， 实 验 结果 如 表 1~3 所 示 。 


environment based on K-means feature [J]. Application Research of 
从 表 1 可 以 看 出 ， 在 道路 区 域 分 割 精度 上 ， 本 文 算 法 比 文 Computers. 2016, 39 (2): 602-606. ) 
献 [7，8] 分 别 高 出 28.1% 和 45.1%， 略 低 于 文献 [9] 的 精度 ; 由 表 [3] 邓 燕 子 , 卢 朝阳 ， 李 静 ,等 . 采用 多 层 图 模型 推理 的 道路 场景 分 割 算法 
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